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Abstract 

We use the canonical coset parametrization and provide a formula 
with the unitary part of the Bures measure for non-degenerate sys- 
tems in terms of the product of even Euclidean balls. This formula is 
shown to be consistent with the sampling of random states through 
the generation of random unitary matrices. 
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1 Introduction 

The Bures measure is the volume element of the corresponding Bures metric 
that can be obtained from the infinitesimal form of the quantum fidelity be- 
tween mixed quantum states [UEIE], or from the statistical distance between 
mixed quantum states [4J. The Bures measure has been proposed as a prior 
distribution for implementation of quantum Bayes estimation [5l E] . Other 
proposals and other measures also exist including monotonic Riemannian 
measures [7J Hj. The term monotonic is applied to stochastic maps, which 



are not allowed to increase the distance. Andai [8] calculated the volume 
of the whole state space according to the Lebesgue measure including a few 
monotonic Riemannian measures [TJ. 

The quadratic form of the Bures metric can be written as 

dB(p,p + dp) 2 = ^Tr[Gdp}, (1) 

with G implicitly defined from dp = Gp + pG. A more practical formula was 
found by Hiibner [9] in terms of the eigenvalues (populations) of the state as 

^(P,P+^) 2 = ^EE l(A f!'y |2 ' 0) 

where p\\k) = Ajt|Afc). A general state can be parametrized by applying a 
unitary operator on the diagonal state p^ = diag(Xi, A2, ...Ajy) as 

p = SlpWtf, (3) 

with the unitary operator [TOl [TT] as the generalized flag manifold 

Qe TT( — n ~ n} N }^ tt? — v m 1 +m 2 + ...m q = N (4) 
(/(mi) <8) U{m 2 ) <g> ...U(m g ) 

where is the degeneracy of the unique eigenvalue Ay. The flag manifold 
can be decomposed as a product of cosets in order to develop a suitable 
parametrization. For example, the non-degenerate case can be decomposed 
as 

o c U(N) _ U(N) U(N-l) U(2) 

U(1)® N U(N-1)®U(1) U(N -2)®C/(1)"'C/(1)®C/(1)' 1 ' 

The eigenstates can be written as |Afc) = Q\k) and the infinitesimal variation 
of p can be expanded as 

dp = Qdp {D) Q j + dVLp {D) & + Qp^dQl (6) 

Introducing this expression in Hiibner's formula (j2j) we find 

4 A, A,- + Afc 



. AT , s 2 V N 

iE^ + EE A iJfc |ntdn|2 fc) (8) 



jr' = l " 3 j=l k=j + l 
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with An- = ^ , , , such that the volume element can be extracted to obtain 

j + k ' 

Hall's formula [12] up to a scale factor. The volume element with the scale 
used in [13J is 

dV B = 5(Ai + A 2 + ...Ajv - IT A jk dx jk dy jkl (9) 

where dxj k = Re(Q Jf dQ)jk and dyj k = Im(Q Ji dQ)j k . The remarkable feature 
of this expression is that it separates the effect of the eigenvalues (popula- 
tions) and the effect of the unitary operator, 

dV B = dV^dVP. (10) 

This paper is organized as follows. Section 2 reviews some results of the Euler 
parametrization and introduces some concepts and formulas to be used later 
in the paper. The main contribution of this paper is developed in section 
3. Section 4 compares the results with the generation of random unitary 
matrices. Section 5 presents concluding remarks. 



2 Euler Parametrization 

The generalized Euler parametrization was developed by Tilma and collabo- 
rators p~U [15] and was used in [16] to calculate the volume and measure of 
the unitary part of the Bures measure. In this section we review the three- 
level case in order to introduce some concepts and formulas that will be used 
in the next section. The unitary operator for a three-level system can be 
parametrized as [151 EH E] 

tt = (e^ A3 e 05A2 e* 4A3 e </,3A5 )(e </,2A3 e </,lA2 ), (11) 

with the factors of f2 as parametrisation of u^^uii) anc ^ u(r)®u(i) m terms 
of the Euler angles <ft k and the Gell-Mann matrices A^. The measure of the 
unitary portion can be calculated as the product of the measure of the corre- 
sponding cosets. For non-degenerate 3-level systems, the coset decomposition 
is 

dV™ = dv( Wl f (3 > Wf„, U(2 l .) (12) 



U(2)®U(1)J \U(l)(g)U(l 
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The volume of the coset u ^2)^u(i) can ^ e obtained by calculating tt^dVl and 
selecting the matrix components where the corresponding Lie algebra lies. 
Thus, we may extract the coordinate transformation from the following terms 

• dx 3 ' + idx 4 '\ 

rtdn = I • • dx 5 ' + idx G ' (13) 



to obtain 



dV ( , , A = cos 3 sin 3 (f) 3 sin 2(f) 5 d<j) 3 d<j) 4 d<j) 5 d<j) 6 . (14) 

\U(2)®U{1)J 

A similar procedure can be carried out for the second coset. 

The Haar measure of the coset EfTn^Tm * s topologically equivalent to an 
even sphere S 2n according to Gilmore [HI [19] with the corresponding volume 



However, the measure of the unitary section of the Bures measure is not 
the Haar measure. Some references refer to it as the truncated Haar Measure 
[15] [20] . Direct integration of the coset measure does not result in the volume 
of even spheres r 2 = 1, but instead in the volume of even balls r 2 < 1, with 
r as the radial coordinate, such that 

(16) 

2W 2 

Vol(iT) = — — -, (17) 
ni (n/2) 

which is in perfect agreement with the formulas found in [211 HB] • 

Consequently, the volume of the unitary section for non- degenerate sys- 
tems is equal to the product of the volume of even balls 



with 



( U(N) 



TV 



N(N-l)/2 



V0 ' [Ufa* ) = Vd(B )Vol(B 2 "-)...Vol(B^) = (18) 
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This result is also consistent with the volume presented by Sommers and 
Zyczkowski [5] as 

, / U(N) \ (2tt) jv ^- 1 )/ 2 



Vo1 = ~wm- (19) 

The discrepancy factor can be explained by a simple numerical scale factor 
of 1/2 on the Bures metric, because N(N — 1) is equal to the dimension of 
the Lie algebra occupied by the coset space. Equivalently, N(N — 1) is the 
number of degrees of freedom required to parametrize jjwjm ■ So far, we have 
shown that the volume of the coset (fT6|) can be written as the volume of an 
even ball, but this does not imply that the measure of the coset is Euclidean 
defined on an even ball. This assertion is proved in the next section and 
further numerical tests are carried out in section 4. 



3 Canonical Coset Parametrization 



An important parametrization arises from the canonical coset, as presented 
by Gilmore [18] on page 351. The Bures metric was obtained for 3-level 
systems in [22J and a more general prescription in [11] for N-level systems, 
but the measure was not calculated in this formulation. The power of the 
canonical coset parametrization lies in the many possibilities to analytically 
express the exponential of the following typical block matrix 



exp 



fi 

-fit 



cos 



VBB^ 



sin vWfB 



B 



VB^B 

sWBtB ^t COS VB^B I 



(20) 

where the case of interest is such that B = column^z 1 , z 2 , ...,z N ~ l ) is a 
column vector of complex numbers. This exponential can also be expressed 
in terms of spherical coordinates x J as 



exp 



B 

-fit o 



[l-XXt] 1 /? 

-xt 



X 

[1-XtX] 1 / 2 



(21) 



such that 



/ 



X 



sin VB^B 



B^B 



-B 



x 1 + ix 2 
x 3 + ix A 



\ 



x %N-3 _|_ j x 2N~2 



(22) 
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This coordinate system is called spherical because the column vector is made 
of variables that range inside an even ball B 2k , where the radial coordinate 
is r 2 = X'X. The exponential in (I2ip is important because it provides a 
parametrization of the coset ^(jv-i)®i7(i) as a N x N matrix. The coset 
required to parametrize the unitary section of the Bures measure can be 
constructed in terms of products of layered cosets flSJ) that have form 



U(2) 



17(1) <g) Z7(l) 



U(3) 



U(2)®U(1) 



IX 



/V 1 - (a: 1 ) 2 - (^ 2 ) 2 __J 

-(x^ix 2 ) y/1 - (x 1 ) 2 - (x 2 ) 2 



-(a; 3 — ix A ) 



w. 



(2) 



r + ix 
x 5 + ix 6 



-(x 5 



; ) y/l - (* 3 ) 2 " 



0.. 


o\ 


0.. 





0.. 





1.. 





0.. 
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0... 


o\ 


0... 





0... 





1... 





0... 


V 



with — — XX^)j k . With this background, we state the following 
theorem 

Theorem 1 The measure of the following coset corresponds to an Euclidean 
measure defined inside of an even ball, such that 



dV 



U{n + 1) 



dV E {B 



2 m 



U(n)®U(l), 
This theorem applied to 3-level systems results in 

17(3) 



dV 



U(2)®U(1) 



dx 3 dx 4 dx 5 dx 6 , 



(23) 



(24) 



such that (x 3 ) 2 + (x 4 ) 2 + (x 5 ) 2 + (x 6 ) 2 < 1, in terms of the variables of the 
corresponding canonical coset parametrization. 

Proof. The strategy is based in the generalization of the proof initially 
provided for the simpler case dV ^ , , . , , , 
interest is 



i ... u( - 3 ).. . ) The complex column X of 



X 



x 3 + ix A 
x 5 + ix 6 ' ' 



(25) 



such that the unitary operator becomes 



n 



1/2 



X 



-Xt [1 - x*x] 



1/2 



i + ^f- x xxt X 



(26) 



with r 2 = X^X = (x 3 ) 2 + (x 4 ) 2 + (x 5 ) 2 + (x 6 ) 2 . The measure is invariant under 
an orthonormal transformation O applied to the coordinates ( 
This means that it is sufficient to consider the evaluation of Q^dQ at ( •JO j tC j OC j OC 1 
(r, 0, 0, 0), which produces the following expression 



( —irdx A 

(y/l^-l)(dx 5 +idx 6 ) 
dx 3 +i(r 2 -l)dx 4 



(VT^-l)(dx 5 -idx 6 ) dx 3 -i(r 2 -l)dx 4 





-dx 5 + idx 6 



dx 5 + idx 6 
irdx 4 



\ 



J 



(27) 



fdx 3 '\ 
dx A ' 
dx 5 ' 
\dx 6 'J 

which leads to the measure 



The coordinate transformation can be extracted from (Q*dQ)i3 and (TvdQ)23 
as 

0\ fdx 3 \ 
VI -r 2 
1 
1/ 









V o 



dx A 
dx 5 
\dx 6 / 



(28) 



dV 



U(3) 



U(2)®U(1) 



dx 3 dx A dx 5 dx 6 , 



(29) 



with r < 1. The transformation matrix (128]) changes with the application 
of an orthonormal operator O on the coordinates ( Ob ■ Ob ■ 0C* ■ 0C* 6 ), but the de- 
terminant remains invariant as stated before. The only differentials without 
trivial transformation are dx 3 and dx A which were evaluated at a; 3 = r and 
x A = 0. The rest of the differentials transform according to the identity. This 
means that by extending the coordinates to higher dimensions the extra dif- 
ferentials will transform according to the identity as well, thus, maintaining 
the determinant equal to 1 and proving the theorem. ■ 

This theorem leads us to formulate the following formula of the Bures 
measure for a state with non-degenerate spectrum 



dV^ n) (N) = dV(B 71 "-' z )dV(B' ZN - i )...dV(B 



->2N-2 



2N-A\ 



(30) 
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Some degenerate states including pure states and those without full-rank 
can be treated by reducing the degrees of freedom and the number of balls in- 
volved in the parametrization. For example, Table [1] shows the characteristic 
diagonal states along with their corresponding measures in low dimensions. 



Diagonal state 


n 


Measure 


diag{0, 0, 1) 
diag{0, 0,0,1) 
diag{0, 0, A 2 , Ai) 
diag(0, 0,0,0, 1) 
diag{0, 0, 0, A2, Ai) 
diag(0, 0, A 3 , A 2 , Ai) 


U(3) 
(7(2)0(7(1) 
1/(4) 

(7(3)0(7(1) 
1/(4) (7(3) 
1/(3)0(7(1) £7(2)0(7(1) 
1/(5) 
(7(4)0(7(1) 
(7(5) (7(4) 
(7(4)0(7(1) £7(3)0(7(1) 
C/(5) (7(4) U(3) 


rf\/(5 6 )rf\/(5 4 ) 
c/y(5 8 ) 

c/V(5 8 )c/V(5 6 ) 
dV(B 8 )dV(B 6 )dV(B 4 ) 


(7(4)0(7(1) (7(3)0t/(l) (7(2)0(7(1) 



Table 1: Unitary part of the Bures measure for degenerate states with 
reduced rank, where Xj 7^ A fc . 



4 Random Sampling 

The results from the previous sections can be used to compare the sampling 
distribution of the Euclidean spheres of the Bures measure against the sam- 
pling distribution obtained from a generation of random unitary matrices. 
The most efficient and transparent method to generate random unitary ma- 
trices is described by Mezzadri |23j, which is based on the QR decomposition 
of complex random matrices. A random state p can be generated by two 
independent methods 

• 1: Through the generation of random unitary matrices. 

• 2: Through the generation of random points on the even Euclidean 
balls B 2k and subsequent use of the canonical coset parametrization to 
obtain the state. 

The alternative distributions seem to be equivalent as can be verified by 
plotting their cumulatives against each other and observing a linear one-to- 
one correspondence up to some fluctuations. A specific test can be designed 
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Figure 1: Plot of the equivalent cumulatives of the component (p)33 for 1000 
random states with spectrum diag (|, |, |). 

for states having the spectrum of the following non-degenerated diagonal 
state 

P W = diag (1,1 ~). (31) 

A plot of the two cumulatives against each other for the (p)33 component is 
shown in Figure [TJ This test was comprehensively carried out and verified 
for systems up to 5-levels. 

5 Conclusions 

We calculated the unitary part of the Bures measure for non-degenerate 
systems in terms of the canonical coset parametrization and found an ex- 
pression as the product of Euclidean even balls. This result was shown to be 
in agreement with the numerical random sampling of unitary matrices and 
with the formulas of the volume found in the literature. These results are 
also relevant to monotone metrics other than the Bures metric including the 
Hilbert-Schmidt measure. 
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